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Abstract 

This  paper  presents  theory  and  analysis  of  single-tip 
field  emission  and  electron  beam  propagation  in  the 
electrostatic  focusing  fields.  It  is  shown  that  two  gate 
apertures  with  a  focusing  anode  allow  transport  of  narrow 
electron  beams  over  long  distances  without  need  for  a 
confining  magnetic  field.  Physical  mechanisms  of  the 
beam  formation,  transport,  field  emission  energy 
distributions,  the  effects  of  the  emission  properties,  and 
parametric  studies  are  discussed,  emission  current 
formula  is  derived,  and  new  concept  and  model  of  the 
bandgap-spread  multilevel  field  emission  is  given. 


I.  INTRDUCTION 

Because  of  their  small  emitting  area  and  high  emission 
current  density,  tip-based  field  emitters  produce  beams  of 
low  intrinsic  emittance  and  high-brightness.  Most 
prospective  applications,  such  as  high-resolution  X-ray 
imaging,  free  electron  lasers,  and  high-frequency  TWT  and 
terahertz  periodic  structures,  require  a  focused  electron 
beam  to  be  transported  over  distances  many  orders  of 
magnitude  larger  than  the  typical  field  emitting  geometry 
transverse  dimensions. 

This  paper  presents  theory,  simulations,  and  analysis  that 
describe  focusing  and  transport  of  the  electron  beam  from  a 
field  emission  tip  in  the  electrostatic  field  produced  by  two 
gate  apertures  and  a  focusing  anode.  The  idea  of  two  gates 
was  considered  first  by  W.  B.  Hermannsfeldt  [1]  and 
further  investigated  by  others  [2-3].  However,  the  two-gate 
concept  does  not  allow  sufficiently  long  focused  electron 
beams.  We  introduced  a  focusing  anode  into  the  geometry 
and  obtained  optimized  configurations  that  provide  the 
desired  long  aspect  ratio  focusing  and  beam  transport.  For 
more  adequate  description  of  the  field  emission,  we  derive 
an  emission  model  with  higher  order  corrections  to  the 
Fowler-Nordheim  theory  and  suggest  a  new  concept  of 
bandgap-spread  multilevel  field  emission.  We  also  discuss 
the  calculated  particles  transverse  energy  distributions, 
consider  effects  of  the  emission  parameters  on  the  beam 
properties,  and  show  an  example  of  a  multiple-beam  field 
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emission  source  that  could  be  integrated  with  a  terahertz 
periodic  structure  for  extended  beam-wave  interaction  with 
no  magnetic  field  for  focusing. 

II.  FIELD  EMISSION 

The  presence  of  electric  field  causes  electrons  to  escape 
from  a  surface  of  the  emitting  matter.  Such  field  emission 
of  electrons  can  occur  into  vacuum  or  a  medium  and  the 
emitting  matter  can  be  electrical  conductor,  semiconductor, 
dielectric,  or  more  complex  substance  in  the  solid  or  other 
state.  In  solid-state  electrical  conductors,  e.g.  metals,  the 
valence  electrons  possess  the  conduction  energy  band  and 
are  described  by  Sommerfeld  free  electron  gas  model  with 
Fermi-Dirac  statistics,  which  defines  the  electrons  energy 
distribution.  For  the  emission  from  not  electrical 
conductors  the  Sommerfeld  theory  of  metals  with  Fermi- 
Dirac  statistics  can  still  be  applied  as  an  approximation  and 
with  suitable  modifications. 

At  sufficiently  high  temperatures  the  energy  of  thermal 
motion  of  the  electrons  allows  them  to  overcome  the 
binding  potential  barrier,  described  by  the  work  function, 
thus  producing  the  thermionic  emission,  which  is  a  special 
case  of  the  field  emission.  The  Maxwell-Boltzmann 
statistics  at  such  high  temperatures  determines  the  electrons 
energy  distribution,  which  is  also  a  special  case  of  the 
Fermi-Dirac  distribution.  The  Richardson-Dushman 
equation  with  the  Schottky  reduction  of  the  effective  work 
function  describes  the  thermionic  emission  current  density. 

The  high  electric  field  triggers  electron  emission 
independently  from  the  temperature.  Such  cold  emission 
occurs  due  to  quantum  tunneling  of  electrons  through  the 
surface  potential  barrier.  The  electric  field  effects  on  the 
shape  of  the  potential  barrier  by  making  it  of  a  finite  width, 
in  the  length  units,  reducing  to  zero  with  the  higher  electric 
field.  Fowler  and  Nordheim  derived  a  theory  of  the  field 
emission  in  1928  [4,5].  Burgess,  Kroemer,  and  Houston  in 
1952  [6]  corrected  a  mathematical  error  in  the  Nordheim 
analysis  of  the  image  charge  effects.  Murphy  and  Good  in 
1956  generalized  the  Fowler-Nordheim  theory  to  describe 
the  cold,  thermionic,  and  intermediate  emission  regimes 
[7].  They  also  suggested  an  improved  form  for  the 
Nordheim  elliptic  function  and  introduced  a  correction  to 
the  Fowler-Nordheim  formula  by  adding  into  the  analysis  a 
linear  term  from  the  Taylor  series  expansion  for  the 
momentum  integral  determined  by  the  Schottky-  Nordheim 
tunneling  potential  barrier  with  the  respect  of  the  image 
charges. 
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The  presented  theory  considers  field  emission  starting 
from  the  energy  distribution  of  the  emitted  electrons  in  both 
normal  and  transverse  directions  to  describe  the  field 
emission  phenomena  including  the  intrinsic  emittance  and 
thermal  effects. 

A.  Normal  and  Transverse  Energy  Distribution 

According  to  Sommerfeld  free  electron  gas  model  the 
number  of  electrons  in  unit  volume  in  the  momentum  range 
dpxdpydpz  is  determined  by  the  Fermi-Dirac  distribution 
function  and  the  total  number  of  the  energy  states  in  the 
corresponding  unit  volume  of  the  momentum  phase  space 

N(j)x,Py,Pz)dpxdpydp^  =  ^^^dpxdpydpz  .  (1) 

1+e  fcr 

Here  p^,  Py,  Pz  are  the  electron  momentum  components  in 
the  Cartesian  coordinates  {X,Y,Z)  respectively, 
X(pr,  P7)  is  in  the  units  of  — --r,  E  is  the  electron 

ry  r z.  j  m3(kg-m/s)3" 

kinetic  energy  measured  relative  to  the  bottom  of  the 
conduction  band,  ju  is  the  Fermi  level,  k  is  Boltzmann’s 
constant,  h  is  Plank’s  constant,  and  the  factor  of  2  arises 
from  the  Pauli  exclusion  principle. 

For  the  electron  velocity  phase  space  {y^,  ^z)  (1)  yields 

N(yx,Vy,Vz)dVxdVydVz  =  dv^dVydVz  ,  (2) 

l+el^ 

where  m  is  the  electron  mass,  N(yx,  Vy,  Vz)  is  in  the  units 
of  /  ,  and  E  =  —(vl  +  vl-\- 

In  cylindrical  coordinates  (x,r,a)  selected  so  that 
Py  =  Vj.  sin  a  and  Vz  =  Vy  cos  a  ,  for  which 
dVydVz  =  VydVyda , 
integration  of  (2)  over  a  from  0  to  In  gives 

N (v„  Vr)dv,dVr  =  4n  Qf  ’  (3) 

1+e  kT 

where  N(Vjr,  Vr)  is  in  the  units  of  — ---  . 

Combining  (3)  with  the  continuity  equation  for  the 
electron  flux  in  the  electron  energy  phase  space  (Ex,E-^) 
find  the  number  of  electrons  striking  a  unit  area  of  the 
emitting  surface,  x  is  the  normal  direction,  per  unit  time 

71-.  71-.  4mn  cLEydEr  ^ 

N  (^E^,  Ej^^v^dE^dE^  —  Ex+Er-ii  ?  (^) 

^  1+e  Ter 


where  NiE^,  E-^)  is  in  the  units  of  and 

E-^  =  mVrl2. 

With  the  account  of  the  transmission  coefficient  T) 
arising  from  the  probability  of  quantum  tunneling  of  the 
electrons  through  the  surface  potential  barrier  obtain  for  the 
emission  current  density  j 

d^jjEx.Er)  _  4nem  Dfe)  .  . 

dExdEr  ~  h'^  ’  V  ) 

^  ^  1+e  kT 

where  e  is  the  elementary  charge,  and  the  designation 

The  transmission  coefficient  T)  is  determined  by  solution  of 
the  Schrodinger  equation.  The  Kemble  solution  [8]  for  the 
parabolic  and  so  the  Schottky-Nordheim  potential  barrier  is 

(6) 

where  Dq  is  the  JWKB  integral 


®o(F;,)  =  e-c  ,  (7) 

with  the  exponent 

G  =j{p  +  (p-  E:^y/^v(y)  ,  (8) 

F  is  the  applied  electric  field,  B  is  the  universal  Nordheim 
constant  given  by  B  =  the  work  function, 

u(y)is  the  Nordheim  function  defined  by  the  integral 

I  .h-p-^dp  ,  (9) 


1  -  p - dp 

p 


p+cp-Ex 


X  ,  andy  is  the  Nordheim  parameter  given  by 


47160  P+(P-Ex 

For  the  strong  potential  barriers,T)o(F;^,)  «  1,  i.e.  G  is  large 
enough,  (6)  approaches  to  a  simpler  form 

.  (11) 
Note  that  with  such  approximation  one  can  let 
T)  =  1  for  Ex  >  p  (Po, 

where  (po  can  be  looked  as  a  work  function  reduced  due  to 
the  Schottky  effect  and  determined  by  Ey  from  (10)  aty=l, 

(Po=(p-  ^Je^F/AnSo  . 

From  (5)  obtain  the  current  density  distributions  over 
electron  kinetic  energy  in  the  normal  and  transverse 
directions,  respectively 


dEr 

Ex+Ei — fd 


djXEr)  _  4nem  j 

dEr  ~  I 

^  /  1+e  kT 

do 

or  upon  integration 

/  rP+4> 

djiEr)  _  4nem  j  I  'r)iEx)dEx  j  r 

W 


+  kT  In 


Here  p  +  0  defines  the  integration  limit,  at  which  the 
transmission  coefficient  T)  ~  1.  In  case  of  (II),  0  =  cp^. 

The  normal  energy  distribution  in  the  form  of  (12,  14)  is 
similar  to  that  obtained  by  Fowler  and  Nordheim  for  the 
cold  emission  calculations.  The  transverse  energy 
distribution  (13,  15)  presented  first  in  our  recent  theory  [9] 
is  key  to  understanding  the  intrinsic  emittance  and  thermal 
effects.  Second  term  in  (15)  shows  the  thermionic  emission 
contribution,  which  becomes  insignificant  for  the  cold 
emission  regimes.  The  following  discussion  is  focused  on 
the  cold  emission  current  density  and  corrections  to 
improve  the  accuracy  for  the  barriers  satisfying  (11). 

B,  Emission  Current  Density 
Integrating  (14)  over  (or  (15)  over  E^  and  using  (II) 
find  the  field  emission  current  density 

rp/(p 


hT\n(l-fe^^d6  ,  (16) 
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where  0((5)  =  (1  +  6y^^v(y)  , 

fi-Ex  _  _  /  _  ^^F/(p 

1+S  ' 


6  = 


(p 


y  = 


and 


(17) 

(18) 


kT\n[l  +  e^^ 


Sq  is  electric  constant.  In  case  of  the  lower  temperature 

5(p 

^\  \  6(p 

IkTekT 

and  (16)  is  simplified  to 

n[l/(p 

4nem  n  I 

^  =  — ^ 


5  >  0,ekT  »  1 

S(p 

S  <0,ei^  «1 


B(p 


3/2 


-0(5) 


8d8 


(19) 


Obtain  Taylor  series  expansion  for  0((5)  defined  by  (17) 

(t)(5)  =  u  +  -T(5  +  -/52  +  0(53)  ,  (20) 

where  v,  t,  and /  are  the  correction  hinctions  designated  for 
brevity  as  v  =  u(yo),  t  =  T(yo),  and  /  =  /(yo), 
respectively,  and  yo  is  the  argument  given  by  (18)  at  ^0, 
yo  =  yj /luEo  ■  yfF / (p  ,  (21) 

(22) 

(23) 


T  =  v--yoV  , 

/  =  n  +  l(yoV'-yon')  , 

,  dv\  1  //  d^v\ 

V  =-\  ,  and  n  =  ^ 

“^'y=yo  ^  y=yo 


Here  the  Nordheim  function  v  is  given  by  (9)  and  can  be 
expressed  through  the  elliptic  integrals  [7] 

»0')  =  VTT7(E(lg)-yK(a)),  (24) 

where  K  and  E  are  the  complete  elliptic  integrals  of  the 
and  2^^  kind,  respectively,  commonly  defined  as  follows 


rn/2 

Jo 


E(t)  =  Jq  ^  Vl  —  t  sin^  8  M, 


K(t)  -  ,  ^ _ _ 

V  l-tsin^ 

|arg(l-t)|  <n . 

Taking  into  account  the  found  derivatives 


dv 

dy 

d^v  _ 

2^fT+y 

3/4 

dy^ 

(l-y)Vl+y 

o 

(E(iS)-o (tS)) 


obtain  from  (23) 

Considerations  similar  in  principle  to  (19-25)  have  been 
conducted  by  Murphy  and  Good  [7].  However  their 
analysis  is  limited  by  the  order  correction  function  r. 
One  can  see  though  from  (17-20)  that  in  case  of  the  smaller 
work  function  g?  (e.g.  reduced  by  the  Schottky  effect  or  the 
band  bending  effects  in  semiconductors,  or  due  to  the 
emitter  properties),  the  variable  S  increases,  which 
degrades  accuracy  of  the  series  representation  (20).  Figure 
1  shows  the  exact  values  of  0((5)  obtained  from  (17)  and 
from  the  expansion  (20)  calculated  for  the  linear,  using  v 
and  r,  and  quadratic,  using  v,  r  and  /  approximations. 
Figure  2  shows  plots  of  the  correction  functions.  Cutler  and 
Good  in  the  follow-up  paper  [10]  have  included  the  f- 
function  into  the  field  emission  calculations,  due  to 
different  motivations.  However,  this  was  done  with 
additional  approximations  by  taking  the  f-tQxm  out  from  the 
exponent  and  other  expansion  terms  to  simplify  the 
emission  current  formula  for  further  integration,  which  in 
turn  reduces  the  accuracy  and  imposes  constraints  on  the 
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Figure.  1.  Exact  values  of  0(5)  from  (17)  versus  its 
linear  and  quadratic  approximation  from  (20).  The  plots 
are  shown  for  the  example  with  F=2  GV/m,  (p=A  eV 
which  correspond  as  well  to  other  electric  field  and  work 
function  values  with  the  same  ratio  of  Vf/ cp. 


Figure.  2.  Correction  functions  v,  r,  and /  The and  '^/f 
dashed  lines  are  shown  for  further  discussion  on  the 
correction  to  the  field  emission  formula. 


allowable  y^values.  We  have  found  that  no  additional 
simplifications  are  necessary  for  the  correct  integration  and 
included  the  higher-order  correction  function /into  the  field 
emission  formula. 

From  (19,  20)  obtain 


j  = 


g3  p2  B(P^F 

- -e  p 

Snh  (pr^ 


■z(d,c) 


(26) 


where  2' (^,  C)  is  the  correction  factor  given  by  the  integral 


Z(ig,C)  = 


e  ^  "^3  vdv 


(27) 


3B  (D^ F 

V  =  — —  T  •  8,  and  g  and  C  are  the  parameters  defined  as 


g 


3B  (p 


3/2  -2 


c  = 


,and5e  =  -^T. 

^  2  F 


(28) 


2  F  f  '  (pT 

The  left  part  of  (26)  without  the  2-factor  is  the  known 
Fowler-Nordheim  formula  with  the  o  (Nordheim  [5])  and  t 
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Figure.  3.  Plots  of  the  function  A(g,  C)given  by  (32). 


Figure.  4.  Plots  of  the  correction  factor  2' (^) determined 
by  (3 1).  The  dots  show  example  points  for  given  F  and  (p. 


(Murphy  and  Good  [7])  correction  functions.  Explicitly,  the 

correction  2-factor  integral  (27)  is 

Z(g,e)  =  2gx  ^(29) 

(l  -  [erfc^  -  erfc  +  f))]  “  . 


It  can  be  simplified  by  presenting  Z(g,C)  in  the  form 
Zig,e)  =  Z(ig)(l-Mg,e))  , 
where  Z (g)  =  2g{l-  erfc ^5) 


and 


A(5<  O 


e  ^^(^'*'4)— erfc^V5(l+f)) 
1-F^ed  erfcyfg 


(30) 

(31) 

(32) 


A  Taylor  series  expansion  of  the  2  (^) -factor  (30), 
originally  found  with  the  symbolic  mathematics  software 
Mathematica,  is  given  by 

2:(^)  =  +  +  ■  (33) 

2g  4a2  8o3  \g^J  ^  J 


Figure  3  shows  plots  of  the  residual  function  A(g,C).  In 
many  cases  typical  value  of  C  is  greater  than  1  and  g  is 
greater  than  3,  determined  by  (28),  so  that  A(g,C)  «  1  and 
the  correction  factor  2 (^)  from  (31)  plotted  in  Fig.  4  can 
be  used  in  the  emission  formula  (26)  instead  of  2(^,  C). 


III.  MODELING  OF  FIELD  EMISSION 

Based  on  the  physical  mechanisms  and  the  number  of  the 
emitting  energy  bands,  in  describing  and  modeling  the  cold 
field  emission  we  distinguish  three  types  of  the  emission. 
These  are  the  single-band  and  the  two-band  field  emission, 
and  the  suggested  new  type  of  bandgap-spread  multilevel 
field  emission. 


A.  Single-Band  Field  Emission 

For  electron  emission  from  a  single  energy  band,  that 
may  occur  for  example  in  electrical  conductors  or 
dielectrics,  we  have  implemented  into  the  MICHEFFE  [11] 
particle  optics  code  the  improved  emission  model  with  the 
correction  to  the  Fowler-Nordheim  theory.  The  correction 
takes  into  account  the  higher  order  series  expansion  terms 
for  the  integral  determined  by  the  Schottky-Nordheim 
tunneling  potential  barrier.  As  discussed  in  the  section 
above  the  approximation  error  without  this  higher  order 
correction  becomes  critical  for  the  lower  work  hinctions 
relative  to  the  energy  spread  of  the  emitted  particles.  The 
current  density  j  emitted  from  each  site  of  the  emitting 
surface  is  described  by  (26)  with  the  correction  factor 
Z{g,  C)  from  (29)  or  its  simplified  form  2(^)from  (31). 

One  can  evaluate  the  Nordheim  function  u,  and  the 
correction  functions  r  and  /  from  (24),  (22),  and  (25), 
respectively.  Alternatively,  to  avoid  evaluation  of  the 
elliptic  integrals  in  (24)  we  can  use  the  series  expansion 
representations.  Dean  and  Forbes  derived  series  expansion 
for  u(0,  where  I  =  [12].  Upon  regrouping  present  their 

series  expansion  in  the  following  ready  to  use  form 


where  ag  =  1  ,  /cq  =  1  +  61n2 


\nl)l 


i+l 


(34) 


i+3/[16(i+l)]  , 

a^+i  = - 77:; - ai  ,  and 


i+2 


fci+i  =  fci  -  7 


i+1/4 


-^  +  — +  — 
1+3/4  i  +  l  1+2 


We  show  here  the  recursive  rather  than  explicit  relations 
for  the  a^+iand  ki+i  coefficients  since  this  is  preferential 
for  implementation  into  the  computer  simulation  process. 
From  (22)  and  (34)  obtain  series  expansion  for  t 

t(0  =  1 +  (fcj -InO  -  ij  •  (35) 


The  correction  function / can  be  evaluated  from  (25)  and 
(24)  or  through  a  series  expansion.  From  (23)  find  that 

f  =  V  +  —  ^  ,  which  in  combination  with  (34) 

3  dl^\l=y^ 


yields 

/(/)  =  ^  (36) 

1  -  ^  Ui  ((i(i  +  1)  +  ^)  (fci  -  InO  -  2i  -  l)  1^+^ . 


B.  Two-Band  Field  Emission 

The  two-band  field  emission  model  suggests  emission 
from  both  conduction  and  valence  energy  band  that  may 
occur  for  example  in  silicon  or  other  semiconductor 
materials.  The  model  applies  the  aforesaid  single-band  field 
emission  model  to  each  of  the  emitting  energy  bands.  The 
two-band  model  takes  into  account  the  field  penetration 
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and  band  bending  effect  on  the  effective  work  hinction  as 
well  as  the  supply  limited  conditions  when  emission,  e.g., 
from  the  conduction  band,  is  saturated  due  to  the  limited 
number  of  electrons,  which  is  controlled  to  a  substantial 
degree  by  doping. 

C.  Bandgap-Spread  Multilevel  Field  Emission 

For  describing  the  field  emission  from  diamond  we 
introduce  a  new  concept  of  bandgap-spread  multilevel  field 
emission,  which  involves  emission  from  the  multiple 
energy  states  that  are  either  discretely  or  continuously 
distributed  within  the  bandgap.  The  bandgap  spread  can  be 
caused  by  various  kinds  of  imperfections,  intentional  or 
accidental,  such  as  lattice  defects,  dislocations,  phase  and 
chemical  impurities,  doping,  and  surface  treatment.  This 
emission  model  can  be  applied  to  other  materials  as  well 
that  emit  electrons  from  multiple  energy  states. 

As  an  example,  the  suggested  bandgap-spread  multilevel 
field  emission  model  explains  the  voltage-current 
characteristics  observed  from  a  single  diamond  field 
emission  nanotip.  We  can  see  in  Fig.  5  good  agreement 
between  the  simulations  using  the  bandgap-spread 
multilevel  field  emission  model  (3 -emission  levels  in  this 
case)  and  experiment.  On  the  contrary,  the  single-  or  two- 
band  models  cannot  explain  these  experimental  results. 


Figure.  5.  Simulations  of  the  electron  gun  with  diamond 
nanotip  using  the  new  bandgap-spread  multilevel  field 
emission  model  (3 -emission  levels)  show  good  consistency 
with  the  experiment.  Experimental  data  shown  are  courtesy 
of  C.  A.  Brau  of  Vanderbilt  University. 


IV.  EFFECTS  OF  THE  FIELD  EMISSION 
PROPERTIES 

It  is  imperative  to  understand  the  effects  of  the  emission 
properties  on  the  electron  beam  focusing  and  attainable 
characteristics.  With  such  understanding  one  can  determine 
and  optimize  requirements  on  the  field  emitters. 
Furthermore,  the  precise  field  emission  properties  may  be 
not  known  or  not  achievable,  so  one  can  describe  the  beam 
focusing  and  propagation  in  the  channel  as  a  function  of 
emitter  parameters.  For  example,  our  simulations  show  that 
the  effects  of  dopant  level  on  emission  properties  tend  to 
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Figure  6.  Effect  of  doping  on  field  emission  current  in  the 
double-gate  geometry  with  focusing  anode  shown  in  Fig.  7. 


saturate  with  larger  dopant  amounts,  as  seen  in  Fig.  6.  Also, 
a  larger  work  fimction  in  general  degrades  the  beam 
quality,  affecting  the  emittance  and  brightness.  However,  in 
most  cases  the  focusing  structure  remains  functional  and 
requires  only  minor  voltage  adjustments  to  accommodate 
the  variation  in  field  emission  properties. 


V.  ELECTROSTATIC  FOCUSING 


The  main  requirements  for  the  focusing  geometry  are 
that  it  should  capture  as  many  as  possible  of  the  emitted 
electrons  and  transport  and  focus  them  with  minimal  loss 
over  the  distance  and  at  the  diameter  required  for  the 
particular  application.  Power  losses  of  the  entire  structure 
should  be  minimized.  In  addition,  the  focusing  approach 
should  be  expandable  to  larger  field  emitter  arrays  and  yet 
remain  feasible  for  fabrication,  taking  into  account  the 
micron  scale  of  the  geometry.  Simulation  results  of  field 
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Figure.  7.  Simulation  of  field  emission  with  electrostatic 
focusing.  Colors  represent  the  current  in  the  beam  annular 
layers.  Top:  Entire  model  with  particle  trajectories.  The 
image  is  stretched  by  xlOO  in  radial  direction.  Middle: 
Overall  view  of  the  focusing  geometry,  in  scale.  Bottom: 
electron  emission  and  beginning  of  the  beam  formation. 
Inset  shows,  magnified,  the  emitting  lO-nm  radius  tip 
including  the  small  but  still  occurring  side  emission. 
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Figure  8.  Electron  beam  current  density  distribution. 


emission  from  a  single  silicon-tip  with  electrostatic 
focusing  are  shown  in  Fig.  7.  We  derived  this  example 
geometry  for  a  10-100  kV  high  resolution  X-ray  imaging 
source.  The  focusing  optics  has  extraction  and  focusing 
gates  positioned  near  the  emitter  tip  and  a  focusing  anode 
downstream  on  the  beam  axis.  The  extraction  gate  1 
controls  the  emitted  current.  The  focusing  gate  2  pushes  the 
expanding  particle  orbits  toward  the  beam  axis.  The 
focusing  anode  bends  the  expanding  particle  orbits,  thus 
allowing  high  aspect  ratio  confinement  and  beam  transport. 
The  current  density  distribution  in  the  anode  plane  is  shown 
in  Fig.  8. 

We  have  discovered  that  the  suggested  double-gate  with 
focusing  anode  geometry  can  form  fundamentally  better 
quality  electron  beam  with  laminar  trajectories.  Such  new 
laminar  beam  formation  regime  requires  the  sharper  tip 
radius,  lower  work  function,  or  shorter  tip  to  anode  length, 
e.g.  1  eV,  2  nm,  and  5  mm,  respectively,  for  which  the 

^  Focusing  anode  Anode  Side  emission 

Aperture  radius  10pm  plane  not  shown 

100  kV  (negligible) 

z=5mm 

Field  emission  nanotip  Axis  Z 
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I  _ _ 1965  V 
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z=0  mm,  [7=0  V 
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Emitted  i  /  j 

current  1  pA  |  j 

1  j 

Gate  1 

Gate  2  |  i 

4.49  V 

2  V  i - - ■! 

Figure.  9.  Simulation  of  field  emission  with  electrostatic 
focusing  with  the  lower  work  function  emitter  of  1  eV. 
Top:  Entire  model  with  particle  trajectories.  The  image  is 
stretched  by  xlOO  in  radial  direction.  Middle:  Overall 
view  of  the  focusing  geometry,  in  scale.  Bottom:  electron 
emission  and  beginning  of  the  beam  formation.  Inset 
shows,  magnified,  the  2-nm  radius  emitting  nanotip. 


simulation  results  are  displayed  in  Fig.  9  versus  the  4.05eV, 
10  nm,  and  9.8  mm  example  shown  in  Fig.  7.  The  sharper 
tip  and  shorter  length  lead  to  higher  electric  field  formed 
over  the  tip  surface  due  to  the  anode,  including  focusing 
anode,  voltage  that  together  with  the  lower  work  function 
reduce  demand  for  the  electric  field  from  the  extraction 
gate  to  achieve  the  required  emission  current.  As  the 
necessary  gate  field,  and  so  the  transverse  component  of  the 
beam  extraction  field,  is  minimized  the  optics  produces 
laminar  electron  beam  of  the  smaller  emittance,  reduced 
spot  size  with  the  sharp  boundary,  and  vanishing  halo. 

Acknowledgements 

We  appreciate  helpful  discussions  with  Kevin  Jensen  of 
NRF  and  Fuis  Velasquez-Garcia  and  Frances  Hill  of  MIT. 


VI.  REFERENCES 

[1]  W.  B.  Hermannsfeldt  et  al.,  “High-Resolution 
Simulation  of  Field  Emission”,  Nucl.  Instrum.  Methods  A, 
Vol.  298,  (no.  39),  pp.  39-44,  (1990). 

[2]  R.  M.  Mobley  and  J.  E.  Boers  “Computer  Simulation  of 
Micro-Triode  Performance”  IEEE  Trans.  Electron  Devices, 
Vol.  38  (no.  10),  pp.  2383-2388,  (Oct.  1991). 

[3]  J.  Itoh,  et  al.,  “Fabrication  of  double-gated  Si  field 
emitter  arrays  for  focused  electron  beam  generation”  J. 
Vac.  Sci.  Technol.  B  13(5),  pp.  1968-1972,  (1995). 

[4]  R.  H.  Fowler  and  F.  W.  Nordheim  “Electron  Emission 
in  Intense  Electric  Fields”  Proc.  R.  Soc.  A.  Fondon,  Ser. 
119,  pp.l73-181,(1928). 

[5]  F.  W.  Nordheim  “The  Effect  of  the  Image  Force  on  the 
Emission  and  Reflexion  of  Electrons  by  Metals”  Proc.  R. 
Soc.  Fondon  A,  Ser.  121,  pp.626-639,  (1928). 

[6]  R.  E.  Burgess,  H.  Kroemer,  J.  M.  Houston  “Corrected 
Values  of  Fowler-Nordheim  Field  Emission  Functions 
v(y)  and^();)”  Phys.  Rev.  Vol.  90,  (no.  4),  pp.515,  (1953). 

[7]  E.  F.  Murphy  and  R.  H.  Good,  Jr.  “Thermionic 
Emission,  Field  Emission,  and  the  Transition  Region” 
Phys.  Rev.  Vol.  102,  (no.  6),  pp.  1464-1473,  (1956). 

[8]  E.  C.  Kemble  “A  Contribution  to  the  Theory  of  the 
B.W.  K.  Method”  Phys.  Rev.  Vol.  48,  (September  15), 
pp.549-561,  (1953). 

[9]  V.  Jabotinski,  J.  Pasour,  K.  Nguyen,  J.  Petillo,  B. 
Eevush,  D.  Abe  “Electrostatic  Focusing  for  a  Field 
Emission  Electron  Source”  Int  Conf.  on  Vacuum 
Electronics  IVEC  2013  Paris,  France,  2pp.,  (May  2013). 

[10]  P.  H.  Cutler  and  R.  H.  Good,  Jr.  “Higher  Order 
Corrections  to  the  Field  Emission  Current  Formula”  Phys. 
Rev.  Vol.  104,  (no.  2),  pp.308,  (1956). 

[11]  J.  Petillo  et  al.  “The  MICHEEEE  Three-Dimensional 
Electron  Gun  and  Collector  Modeling  Tool:  Theory  and 
Design”  IEEE  Trans.  Plasma  Sci.  Vol.  30,  (no.  3),  pp. 
1238-1264,  (June  2002). 

[12]  R.  G.  Forbes  and  J.  H.  B.  Deane  “Reformulation  of  the 
standard  theory  of  Fowler-Nordheim  tunneling  and  cold 
field  electron  emission”  Proc.  R.  Soc.  A,  Ser.  463,  pp. 
2907-2927,  (2007). 


1006 


6 


